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I.  INTRODUCTION 


The  purpose  of  this  pape7-  is  to  summarize  work  of  Karr  (1972a, 

1972b,  1973)  concerning  a  class  of  models  of  combat  attrition 
processes  which,  in  the  sense  that  attrition  occurs  instantaneously, 
do  not  evolve  over  time  (a.lthough  discrete  time  dynamic  processes 
can  be  constructed  from  these  essentially  static  models  as  we  briefly 
discuss).  We  refer  the  reader  to  Karr  (1974)  for  a  detailed  dis¬ 
cussion  of  a  class  of  continuous  time  parameter  stochastic  processes 
analogous  to  F.  W.  Lanchester's  differential  equation  models. 

Implementation  of  these  static  models  in  computerized  combat  simu¬ 
lations  is  sometime by  means  of  an  exponential  approximation  to  an 
equation  for  expected  attrition  which  is  binomial  in  form  (hence  the 
title  of  the  paper).  We  discuss  the  desirability  and  applicability  of 
such  approximations,  but  only  after  first  giving  rigorous  derivations, 
from  carefully  and  precisely  stated  probabilistic  assumptions,  of  the 
expected  attritions  in  homogeneous  and  heterogeneous  cases  of  point 
fire  combat  (where  a  target  must  be  detected  before  it  can  be  attacked) 
and  for  area  fire  directed  at  a  homogeneous  set  of  targets.  The  meaning 
and  appropriateness  of  the  families  of  underlying  assumptions  may 
then  be  discussed  in  a  useful  manner.  In  addition  to  exponential 
approximations  we  consider  seme  computational  simplifications  in  the 
heterogeneous  point  fire  case  and  difficulties  which  occur  when,  as 
happens  in  an  iterative  computerized  simulation,  random  variables 
are  replaced  in  computations  by  their  expectations. 

The  attrition  equations  we  derive  here  are  one-sided  in  the  sense 
of  expressing  attrition  to  each  force  in  terms  of  its  initial  strength, 
the  initial  strength  of  the  opposition,  and  parameters  describing  the 
physical  situation.  All  attrition  is  suffered  instantaneously.  Two- 
sided  models  can  be  constructed  but  only  in  a  rather  contrived  manner. 
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II.  HOMOGENEOUS  POINT  FIRE  MODEL 


1.  Basic  Derivation 


Consider  a  one-siaed  combat  between  two  homogeneous  forces,  a 
force  of  R  indistinguishable  "targets”  and  a  force  of  B  indis¬ 
tinguishable  "searchers".  We  make  the  following  assumptions  coi  - 
cerning  this  combat. 

Al)  At  a  fixed  time  all  R  targets  become  vulnerable  to 
detection  and  attack  by  the  3  searchers; 

A2)  The  probabiiity  that  the  ith  searcher  detects  the 
target  is  d  for  all  i  =  1,  . . . ,  B  and  j  =  1,  . . . ,  R.  Each  par¬ 
ticular  searcher  detects  different  targets;  independently  of  one 
another; 

A3)  A  searcher  who  makes  no  detections  makes  no  attack.  A 
searcher  who  makes  one  or  more  detections  chooses  one  target  to 
attack  according  to  a  uniform  distribution  over  the  set  of  targets 
he  has  detected,  independent  of  his  detection  process; 

A4)  The  conditional  probability  that  a  searcher  kills  a  target 
given  detection  and  attack  is  k  ,  for  all  searchers  and  targets; 

A5)  No  searcher  may  attack  more  than  one  target; 

A6)  Detection  and  attack  processes  of  different  seai’chers  are 
mutually  independent. 

We  begin  by  computing  the  expected  number  of  targets  killed. 

(1)  PROPOSITION .  Assume  Al)  -  A6)  are  satisfied.  Then  if  K  denotes 
the  number  of  targets  killed,  we  have 

E[K]  =  R{1  -  (1  -  |  Cl  -  (1  -  d)R])B}  . 
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PROOF.  By  elementary  conditional  probability  arguments, 

P{ searcher  i  kills  target  j  } 

=  p{  searcher  i  detects  and  attacks  and  kills  target  j} 

=  P  |  searcher  i  kills  target  j  | searcher  i  detects  and  attacks 
target  j} .  Pjsearcher  i  detects  and  attacks  target  j} 

=  k  P  {  searcher  i  detects  and  attacks  target  j  } 

=  k  P  {  searcher  i  attacks  target  j  |  searcher  i  detects  target  j} 

•  P{ searcher  i  detects  target  j [ 

=  kd  p|  searcher  i  attacks  target  j  |  searcher  i  detects  target  j(  . 


By  the  Law  of  Total  Probability, 

p{  searcher  i  attacks  target  j  J  searcher  i  detects  target  j} 

R-l 

=  v  P  |  searcher  i  attacks  target  j  and  detects  m  other  targets | 

!'  searcher  i  detects  target  j} 

R-l 

=  £  P  {  searcher  i-  attacks  target  j  I  searcher  i  detects 

m-^  target  j  and  exactly  m  other  targets} 

•  P  {  searcher  i  detects  exactly  m  other  targets | 

searcher  i  detects  target  j } 

by  the  identity  P{A  n  B|c|  =  pJa(b  n  C}P|b|c|  . 

By  assumption  A3) 


P{ searcher  i  attacks  target  j  j  searcher  i  detects  target  j 
and  exactly  m  other  targets}  =  —  —  v 

while  by  A2) 

=  P {searcher  i  detects  exactly  m  of  the  other  R-l  targets} 

_  /r-1\  -in  jrR-l-m 
-  )  d  (1  '  d) 
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Therefore 


P | searcher  i  attacks  target  j  j  searcher  i  detects  target  j  j 


=  *£  -i-  l*'1)  d1 

Jo  m  +  1  V  ">  / 


»rn  *  i  j  ■«. R- m-1 

(1  -  d) 


=  V  1  (R-Di  dm 

m  +  1  mi(R-l-m)! 


(1  -  d) 


R-m- 


1 


1  R-  ^  R 1 

*  cR  Jj0  (5  +  1)HR  -  (m  +  1))! 


dm+1(l  -  d)R-(”+1) 


=  A  [i  -  a  -  d>E]  . 


We  thus  obtain 

k  R 

P | searcher  i  kills  target  j}  =  ^[l“vl“d)l 

which,  as  one  expects,  is  independent  of  i  and  j  .  It  follows  by 
the  independence  assumption  A6)  that 

(2)  P | target  j  is  killed} 

=  1  - jp  target  j  is  not  killed} 

B  .  .ii 

-  1  .  p|  fl  jsearcher  i  does  not  kill  target  g  j } 

i=l 

B  . 

—  i  _  n  {P  searcher  i  does  not  kil-l  target  gf 

i=l 

=  i  -  n  a  - 1  [i  -  a  -  d>EU 
•  _  *1  ^ 

1=1 

=  i  -  a  -  |  [i  -  a  -  d)R])B 
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so  that  we  finally  have 


U 

EDC]  =  S  P {target  j  is  killed} 
j=l 

p 

=  E  {l  -  (1  -  |  [1  -  (1  -  d)R]  )B } 
j=l  K 

=  R}i  -  (  i  - 1  ri  -  (i  -  d)Ri)B} 


as  asserted.  f[ 

The  possibility  that  two  or  more  searchers  "kill”  the  same 
target  is  not  excluded. 

2 .  Probability  Distribution  of  the  Number  of  Kills 

To  compute  the  probability  distribution  of  the  number  of  targets 
killed,  which  is  needed  to  extend  this  static  model  to  a  discrete  time 
dynamic  model,  we  take  an  alternative  approach.  It  is  not  true,  as 
one  might  naively  conjecture,  that  the  number  of  targets  killed  is 
binomially  distributed  with  parameters  R  and  q  ,  where  q  is  the 
probability  given  on  the  right  hand  side  ol  (2).  The  reason  for  this 
is  that  even  though  searchers  onerate  ‘ independently  of  one  another, 
different  targets  do  not  die  independently  of  one  another.  Since  each 
searcher  can  attack  at  most  one  target,  knowing  that  a  particular 
target  was  killed  means  that  some  other  target  was  not  attacked  by 
the  searcher  that  killed  the  former  target  and  is  hence  less  likely 
to  have  been  killed. 

(3)  PROPOSITION.  For  m  such  that  m  <  R  and  m'<  B  we  have 

c-ir^ra-qR)  + 

qR  -  kri  -  (i  -  d)R]  . 
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(4)  p|k  =  m}  =  /r 


.  .  si; 

/r\  r 

yn J  r=0 


where 


PROOF.  To  begin,  the  Law  of  Total  Probability  implies  that 
B 

P{K  =  m}  =  2  p{k  =  m,  &  "fatal”  attacks  are  made} 

£=m  ! 


B 

=  E 
£=m 


PjX  =  m|jfc  "fatal"  attacks  are  made} 

P | SL  "fatal"  attacks  are  made} 


where  the  sum  is  zero  if  the  lower  limit  exceeds  the  upper.  Here  a 
"fatal"  attack  is  one  which  would  kill  the  target  attacked  if  no 
other  attacker  attacked  the  same  target.  The  fact  that  two  or  more 
searchers  may  simultaneously  attack  and  "kill"  the  same  target  means 
that  there  may  be  more  "fatal"  attacks  than  targets  killed. 

Now  by  (2),  since  each  searcher  can  kill  at  most  one  target, 

p| searcher  i  makes  a  "fatal"  attack} 

R 

=  E  p{ searcher  i  kills  target  j} 

j=l 


(5) 


=  -  |  [1  -  (1  -  d)R] 
j-- 

=  kf  -  (l  -  d)R] 


But  by  the  independence  assumption  A6)  different  searchers  make 
"fatal"  attacks  independently  of  one  another,  so 


pj exactly 


SL  "fatal"  attacks  are  made} 
(*)  • 


Next,  one  notes  that  computing  the  conditional  probability 
pj  K  =  m| exactly  SL  "fatal"  attacks  are  made}  is  equivalent  to  the 
problem  in  combinatorial  probability  of  finding  the  probability 
that  if  Si  indistinguishable  balls  are  placed  into  R  indistin¬ 
guishable  boxes,  independently  and  according  to  a  uniform 
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distribution,  exactly  R-m  boxss  remain  empty.  Tb.s  is  so  in  view 
of  the  fact  that 

Pj  searcher  i  kills  target  j|  searcher  i  makes  a  "fatal'*  attack}  =  - 

'  R 


for  all  i  and  j  ,  which  follows  at  once  from  expressions  above.  The 
occupancy  probability  is  well-known  and  by  Feller  (1967,  p.  60)  is 


j,  -  »*  o 


(m  -  v)J 


The  Proposition  now  follows  by  some  elementary  computations.^ 

We  may  use  Proposition  (3)  to  construct  a  two-sided,  discrete¬ 
time,  dynamic  version  of  the  unilateral,  static  attrition  process.  ' 
studied  in  Proposition  (1).  We  begin  by  noting  that  we  can  use 
hypotheses  Al)  -  A6)  to  compute  attrition  to  the  B  searchers  if 
roles  ware  reversed  and  these  were  subject  to  detection,  attack, 
and  destruction  by  the  R  targets.  Detection  and  kill  probabili¬ 
ties  are  in  general  a  function  of  the  side  attacking,  so  if  d',  k' 
are  the  parameters  of  the  complementary  process  and  L  denotes  the 
number  of  searchers  killed,  then  by  Proposition^ ?>) 

(6)  P|L  =  n|  =  jo  C-l)n-^)rU-qB)  ♦ 


where 


qB  =  k'[l  -  (1  -  d')B]  . 

One  way  to  realize  a  two-sided  attrition  process  from  our  one¬ 
sided  model  is  to  suppose  that  each  side  is  vulnerable  to  the 
opposition’s  initial  strength,  and  that  the  two  siaes  independently 
attrit  each  other  according  to  equations  (4)  and  (6).  This  is 
unsatisfying  from  a  physical  standpoint  because  it  essentially 
requires  assuming  that  all  detections  and  target  choices  are  made 
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and  then  all  shots  fired  simultaneously. 

This  assumption  does,  however,  lead  to  a  dynamic  model  based 
on  the  stated  assumptions.  The  model  is  a  discrete  time  stochastic 
’-rocess  constructed  on  the  assumption  that  at  each  time  point  an 
interaction  occurs  as  described  in  the  preceding  paragraph,  in¬ 
volving  only  those  searchers  and  targets  that  have  survived 
previous  interactions  and  which  is  otherwise  independent  of  the 
prior  history  of  the  process.  Let  P-.(Ij,R;  k)  be  the  probability 
that  exactly  k  targets  survive  a  single  interaction  if  there  are 
B  searchers  and  R  targets  present;  that  is  P1(B,R;  k)  is 
P|K  =  R-k[  as  given  by  (4).  Similarly,  let  P2(R,B;  j)  be  p{l  =  B-j} 
as  given  by  (6).  Let  Tn  denote  the  number  of  target?  surviving 
after  n  interactions  and  Sn  the  corresponding  number  of 
surviving  searchers. 

The  preceding  statements  together  imply  the  following  result. 

(7)  THEOREM.  Under  the  preceding  hypotheses,  the  process  ((SnJT^))n>1 

is  a  two-dimensional  Markov  process  with  state  space  {o,l,2,„..  \  * 

{ 0,1, 2,... |  and  transition  matrix  P  given  by 

P((B,R),  (j,k))  =  P1(B,R;  j}  •  P2CB,R;  k)  . 


5.  Exponential  Approximations 
Use  of  the  approximation 

C8a)  E[X]  ~  R[1  -  exp(  -  ^[1  -  (1  -  d)RJ)] 

or  the  further  approximation 

-dR 

r  8b)  E[X]  ~  R[1  -  exp(  “  ^C1  “  e  3)3 


Ls  valid  in  the  sense  that 


lim  |e"an  -  (1  -  a)n|  =  0 
n-*»,  a  10 


but  largely  unnecessary,  especially  within  the  context  of  computerized 


combat  simulations.  The  correct  expression  given  in  (1)  can  be 
computed  as  quickly  (or  perhaps  even  more  quickly,  depending  on 
how  a  specific  computer  performs  exponentiations)  as  either  (8a) 
or  (8b).  Moreover,  the  approximations  may  be  rather  poor  for  small 
values  of  B  or  R  ,  or  moderately  large  values  of  k  and  d. 

4.  Use  of  Expectations  as  Inputs 

Another  problem  involved  in  the  use  of  this  model  in  com¬ 
puterized  simulations  is  the  replacement  of  random  variables  in 
certain  expressions  by  the  expectations  of  those  random  variables. 

To  be  more  specific,  suppose  in  tht  Markov  process  of  Theorem  (7) 
we  wanted  to  compute  E[T2],  the  expected  number  of  targets  surviving 
two  interactions.  Let  us  assume  for  simplicity  that  only  targets 
can  be  killed.  Then  subject  to  the  initial  conditions  of  B 
searchers  and  R  targets  we  have 

[:%]  =  za  P2((B,R);  (j.v.)) 

2  j  k 

n 

where  P  is  the  square  of  the  matrix  P  defined  in  (7).  One  might 
seek,  however,  especially  in  simulation  models,  to  use  instead  the 
" approximation” 

e[t2]  ~  r  -  ecTj] 

-  ECTpjl  -  (1  -  -gp^j  [i  -  U  -  d)E[TlV|  , 

with  PTTj]  computed  using  Propostion  (1).  What  is  done  here  is  co 
suppose  there  are  exactly  EfT^  targets  surviving  after  the  first 
interaction,  even  though  this  number  need  not  be  an  integer  (which 
doesn't  matter  too  much,  since  the  right-hand  side  of  (1)  makes 
sense  for  any  B  and  R)  and  to  assess  the  expected  number  of  kills 
in  the  second  interaction  using  E  [T -^]  as  the  initial  number  of 
targets . 
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This  technique  is  simply  wrong;  it  is  not  an  approximation 
because  we  cannot  estimate  the  error  committed  by  its  use,  nor 
give  a  limiting  process  leading  to  no  error,  even  though,  as  the 
example  below  indicates,  that  error  may  be  rather  small. 

The  following  example  gives  an  illustration. 


EXAMPLE,  Suppose  that  B  =  1.  Then  (T  )  is  a  Markov  process 
with  transition  matrix  P  given  by 


P(R,j) 


ik[l  -  (1  -  d)R] 

1  -  k[l  -  (1  -  d)R] 


j  =  R  -  1 
j  =  R 
otherwise 


Suppose 

d  =  .1 
k  =  .5 
Tq  =  3 

Restricted  to  {0,1, 2, 3}  P  is  given  by 


1  0 

.05  .95 

0  .095 

0  0 


0  0 

0  0 

.915  0 

.1405  .8595 


By  direct  computation 


P 


2 


1  0 

.10  .90 

.005  .185 

0  .014 


0  0 

0  0 

.81  0 

.246  .74 


Thus 

ECT-j]  =  2P(3, 2)  +  3P(3,3)  =  2.85 
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and 


ErT2]  =  P2(3,l)  +  2P2( 3, 2)  +  3P2(3,3)  =  2.72  . 

If  we  try  to  compute  E[T2]  as 

ECTild  -  Ej^-j  [1  -  Cl  - 

we  obtain  2.70  which  is  near  to  the  correct  value  2.72  in  absolute 
terms  and  roughly  one  percent  different  in  relative  terms. 

Hence  in  this  case  the  error  is  small,  but  the  error  grows 
with  n  and  cannot,  so  far  as  we  have  been  able  to  investigate, 
be  estimated  or  bounded  a  priori.  Thus  there  remains  the  possibility 
of  substantial  error  if  (1)  is  used  as  an  assessment  equation  in  a 
dynamic  combat  model. 
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III.  HETEROGENEOUS  POINT  FIRE  MOREL 


1.  The  Basic  Model 

In  this  section  we  extend  the  model  of  Section  II  to  the  case 
where  there  are  several  types  of  targets  and  searchers,  with  detec¬ 
tion  and  kill  probabilities  dependent  on  the  type  of  target  and  type 
of  searcher.  The  physics  of  the  process,  however,  remains  unchanged. 

Let  us  assume  that  there  are  M  types  of  searchers,  Bi  searcher 
of  type  i  (i  =  1,  M),  N  types  of  targets,  and  R^  targets  of 

type  j  (j  =  1,  . ..,  N).  We  will  impose  the  following  hypotheses: 

Al)  At  a  fixed  time  all  targets  become  vulnerable  to  detection 
and  attack; 

A2)  The  probability  that  a  given,  fixed  searcher  of  type  i 
detects  a  given,  fixed  target  of  type  j  is  d... 

A3)  Of  the  targets  (of  alx  kinds)  detected  by  a  given  searcher, 
he  chooses  one  to  fire  upon  according  to  a  uniform  distribution. 

A4)  Given  that  he  detects  and  fires  upon  a  target  of  type  j, 
a  given  searcher  of  type  i  destroys  that  target  with  probability  k^; 

A5)  A  given  searcher  detects  different  targets  independently  of 
one  another; 

A6)  No  searcher  may  fire  more  than  once; 

A7)  The  detection  and  firing  processes  of  all  the  searchers  are 
mutually  independent. 

We  wish  to  compute,  under  these  assumptions,  the  expected  number 
of  targets  of  each  type  destroyed.  First,  we  give  an  analytical 
solution.  Let 

B  =  B,  +  ...  +  Bw 
x  M 
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Preceding  page  blank 


denote  the  total  number  of  searchers  and 


R  —  .  +  R^ 

the  total  number  of  targets.  Let  }E  be  the  number  of  type  j  targets 
destroyed. 

(9)  THEOREM.  Under  the  assumptions  A 1)  -  A7)  the  expected 
number  of  type  j  targets  destroyed  is  given  by 

/  M  Bi\ 

(10)  E[K3]  =  -  n  [1  -  dyk^PCiJ)]  J 


where 


K,  ,  R.-l  R. 
j-1  3  0+1 


pd.3)  =  E  ...  E  EE 


ri=0  rj-r°  rr°  rj+r°  rN=0 


- ir — M<&i 

1  +  E  r  \  r .  / 


R.-l-r . 

d,,)  J  J 


P  3 


n  fRqVi?(i 

i<q<N  \r  /  q 


-  d. _) 


R  -r 

^  q  q 


PROOF.  To  begin,  let  us  define  the  following  events 

K^(k,ja)  -  |  target  of  type  j  is  killed  by  kth  searcher  of 
type  i} 

A  -  .  (k,j&)  -  j  jS*"*1  target  of  type  j  is  attacked  by  kth  searcher 
of  type  i} 

D.  .(k,jO  =  { target  or  type  j  is  Etected  by  k  searcher 
of  type  i} 

for  i  <  M,  j  <  N,  k  <  Bi  and  i  <  R . . 
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By  A2)  and  A4),  P|^^(k,£)|  is  independent  of  k  and  £ 
Moreover,  A7)  implies  that 

target  of  type  j  is  killed} 

M  ,  B. 

=  1  -  n  [1  -  P{K ..(1,1)}]  1  . 

i=l  1J 

We  now  need  to  find  the  probabilities  P|K.  .(k,£)[.  To  begin, 

since 


Kjj(k, A)  c  Ai;j(k,jO  c  Di;jtk,je) 

we  have 


HKy(V,*)|  =  P)Kij(k,i)|»ij(k,/)|  •  pjAyCk.jDl 
*  kij  •  pl*ljO‘,i)| 

“  kij  p|Iiij(k.»ll>ij(k,2)|P|Di:jCk>i)| 
=  ky  P|Ay(k,jO|Bij(k,jl)|dij  . 


We  may  write  the  remaining  conditional  probability,  which  is 
independent  of  k  and  &  for  fixed  i  and  j  ,  as 

PlAyCk.jOlDyCMX 


V1/  El 

E  E  ... 

r=o\ri=C 


R.  ,  R.  . 
3-1  3+1 

y  y 

r .  =0  r. 

3-1  3J'J- 


p| j (k,X) (kj£)  n  G^j(r^,  •••>  k,£){ 

x  pjQy (*?]_,  ••  •»  v  k>*)i  » 

where  6..(r_,  ...,  r..;  k,£)  is  the  event  nthe  k^  searcher  of  type 
i  detects,  in  addition  to  'che  £  target  of  type  j,  r.  other  type  j 
targets  and  targets  of  type  p  £  j,”  whose  probability  is  independent 
of  k  and  £  .  By  A3) 
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lN5 


1  + 


1 

N 

E 

P=1 


r 

P 


> 


while  by  A5) 


^®ij^rl*  ***5  rjij’  R>^)} 


R. 

d.  .)  3 
il ' 


-1-r . 
3 


x 


(1 


diq 


R  -r 

)  q  q 


y 


which,  as  expected,  is  independent  of  k  and  z  when  i  and  j 
are  fixed. 

Collecting  terms  now  yields  the  asserted  result.  [[ 

Instead  of  the  assumptions  Al)  and  A3)  concerning  the  physics 
of  the  process  of  detection  and  attack,  we  could  assume  that  targets 
become  vulnerable  to  a  given  searcher  sequentially  in  a  randomly 
chosen  order,  with  all  R!  orders  equally  probable,  that  each  target 
is  either  detected  or  not  and  that  the  first  detected  target  is  fired 
upon.  This  process  occurs  once  for  each  searcher  independently  of  the 
processes  corresponding  to  other  searchers,  but  no  targets  are  destroyed 
until  the  end  of  the  entire  process.  This  interpretation  is  useful  in 
considering  allocation  of  fire  problems. 


2.  Fire  Allocation 

It  should  be  noted  that  this  model  has  no  provision  for  searcher 
determination  of  whether  a  given  detected  target  is  to  be  attacked  or 
not.  In  reality  a  searcher  might  be  tempted  to  pass  over  low  value 
targets  in  the  hope  of  detecting  a  high  value  target.  Or,  if  all 
targets  are  vulnerable  simultaneously,  a  given  searcher  would  not 
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choose  among  them  uniformly,  but  would  instead  choose  the  target  whose 
destruction  entails  the  largest  reward  to  him. 


While  the  model  makes  no  provision  for  such  a  choice  mechanism, 
it  could  be  incorporated  in  the  "simultaneously  vulnerable"  inter¬ 
pretation.  Referring  to  the  proof  of  (  9  ),  we  see  that  the  uniform 
target  choice  is  manifested  solely  in  the  conditional  probability 


P)Aij(k,U|Gij(r1. 


1  + 


1 

N 

I 

P=1 


Hence  a  different  rule  for  target  choice  can  be  incorporated 
simply  by  changing  the  form  of  this  conditional  probability,  subject 
to  some  obvious  regularity  conditions, 

EXAMPLE .  Suppose  the  value  to  a  type  i  searcher  of  destroying 
a  type  j  target  is  u . . .  Then 


=  k . .  u .  . 
13  13 


is  the  expected  return  from  an  attack  upon  type  j  target.  The  rule 
"Of  the  targets  detected  fire  at  one  whose  destruction  is  of  maximal 
value"  leads  to 

....  V  k,i)| 


■= -  if  v.  .  >  v.  for  all  m  ^  j  such  that  r  >  0  , 

1  +  Tj  —  lm  J  m  ’ 

0  otherwise. 

Here  we  have  assumed  for  simplicity  that  v^  i-  v^m  whenever  j  /  in. 
Other  examples  can  be  described  similarly. 

3.  Further  Computations 

To  implement  (10)  as  the  attrition  equation  in  a  computerized 
combat  model  would  be  laborious,  and  lead  to  a  cumbersome  result. 
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Let  us  consider  the  difficulties  more  carefully  beginning  with 
some  simple  cases.  First  suppose  that 

N  =  3 


Rl  ~  R2  =  =  i 


In  this  case,  (11)  yields 

1  1 


=  /o  E0  TTTTS  [(1)  di2 

(i)di3  O  -  di3>1"“] 


=  l[d  -  di2)(l  -  d.3)l  • 

+  7F(1  "  d,-o)d-^  +  d-> 9  C1  "  di03 


+  2'^"ui2;  i3T  i3' J 

+  I  Cdi2  di3^ 

=  l-i(di2  +  di3)+l(di2di3)  > 

which  doesn’t  look  too  ugly;  analogous  expressions  for  p(i,2)  and 
p(i,3)  clearly  exist.  If 

N  =  4 


we  have 


R1  =  R2  =  R3  =  R4  =  1 


p(i,l)  =  1  -  \  Cdi2  +  di3  +  di4> 


+  1  ^di2  di3  +  di2  di4  +  di3  di4^ 
"  k  di2  di3  di4  * 


Thu  general  pattern  is  by  now  evident;  we  summarize  it  in  the 
following  result,  whose  inductive  proof  we  omit. 

(12)  PROPOSITION.  For  any  N  ,  if 

R1  =  R2  =  R3  =  ...  =  =  1  , 
then  for  1  <  i  <  M  and  1  <  j  <  N, 


P(i,j) 


- 1 ,  Siif  I  I 

m=2  m  i1  i2>i1 

i2^j 


L  d.  .  d.  .  ...d.  . 

i  >i  1  1,1i  1 }12  1 >1n 
m  m-1  ^  n 


Proposition  (12)  is  not  so  useless  as  it  first  appears:  in  a 
computer  implementation  one  could  assume  that  there  is  only  one 
target  of  each  typo — even  though  seme  targets  have  the  same  prob- 
bilities  of  detection  and  kill  given  detection — then  carry  out  the 
attrition  calculation  using  (10)  and  (12).  If  targets  were  grouped 
into  classes  of  indistinguishable  (except  by  artificial  notation) 
objects,  attrition  AR^  to  objects  in  the  kth  class  would  then  be 
given  by 


ARk  = 


/  M 

i  -  n 

\  1=1 


where  C,  is  the  set  of  indices 
th 

to  the  k  class.  For  fixed  i 
dent  of  j  e  C^. 


til 

j  such  that  the  j  object  belongs 
and  k  ,  d.  .  and  k.  .  are  indepen- 


While  (13)  could  be  implemented  on  a  computer,  its  use  might 
involve  long  and  time-consuming  computations,  especially  if  R  were 
of  the  magnitude  of  the  number  of  soldiers  in  a  moderate- sized  land 
battle.  Hence  it  is  worthwhile  to  seek  simplified  versions  of  (10), 
several  of  which  we  now  proceed  to  consider. 
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The  following  result,  which  aids  in  developing  simplified  versions 
of  (10),  shows  that  the  complicated  form  of  (10)  arises  from  having 
more  than  one  type  of  target,  rather  than  from  having  more  than  one 
type  of  searcher. 


(14)  PROPOSITION.  If  N  =  1,  then 

M  /  k. 


B. 


(M  /  K.  .  1, 

i  -  n^i  -  -2  i>  -  (i  -  d.)R]j  j 

where  subscripts  i  denote  the  type  of  searcher. 
PROOF.  From  (11)  we  have  for  each  i 

R-l 


vR-l-r 


P(i)  =  P(i,l)  =  r^(r  Jd^'di)1 

which,  as  indicated  in  the  proof  of  Proposition  (1),  is  equal  to 

htr(i  "  (1  ~  di)R)  » 


yielding  (15)„  [J 

We  next  discuss  three  simplified  versions  of  (10)  for  the  case 
N  >  1.  The  notation  AR^  =  E[Kj]  is  used  hereafter. 


a.  Successive  Application  of  All  Searchers  to  Each  Type  of 
Target .  If  there  were  only  type  j  targets  then  by  (15)  the  expected 
number  of  such  targets  destroyed  would  be 


(16) 


r  M  / 

i  _  n  I 

t  k.  .  r 

1  -  -21 

L  i=r 

R:[ 

(1  - 


One  means  of  attrition  assessment  would  calculate  each  ARjby  (16), 
but  this  method  clearly  gives  an  advantage  to  searchers  which  is 
unwarranted  in  terms  of  assumption  Al).  In  effect,  this  method 
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allows  aach  searcher  one  shot  at  each  type  of  target,  rather  than 
one  ^hot  altogether. 

b-  Calculation  With  Weighted  Probabilities.  For  each  i  ,  the 
quantity 

-  1  N 
di  ’  5  &  <1* 

represents  a  detection  probability  for  t'’pe  i  searchers  which  is 
averaged  with  respect  to  the  numbers  of  targets  present,  while 


ki  =  5 


is  a  similarly  averaged  conditional  probability  of  kill  given  detection, 

We  could  then  approximate  the  original  attrition  process  by  one 
with  R  targets  of  a  single  type  and  for  each  i  Bi  searchers  of 
type  i  with  detection  and  kill  probabilities  cL  and  k^,  respectively. 
Using  (15),  attrition  to  these  R  targets  would  be 

ar = 4  -  j  4  -  xt  -  -  3i>R]) l)  • 

Against  type  j  targets  would  then  be  assessed  the  fraction  R^/R  of 
this  attrition  AR  ,  so  that 

(17)  AR.  =  R.(l  -  |  [l  -  (1  -  H.rj)  • 

In  (17)  it  was  unnecessary  to  use  the  averaged  conditional 
probabilities  of  kill  given  detection  and  attack,  and  we  could 
as  well  (and  more  accurately)  have  written 

(17-)  AKj  -  s^l  -  ,nx  (l  -  ^A[i  -  (1  -  <*,)*]  )  j. . 

It  is  possible,  using  one  further  averaging  step,  to  use  the 
homogeneous  equation  directly.  The  numbers 
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T 


and 


d  = 


1 

B 


M 

v 

i=l 


h*± 


k  = 


M 

S  B.K. 
i=l 


i  1 


are,  respectively,  probabilities  of  detection  and  kill  given  detection 
averaged  with  respect  to  both  targets  and  searchers. 

We  may  approximate  the  original  attrition  process  by  an  attrition 
process  with  R  targets  of  a  single  type,  B  searchers  of  a  single 
type,  ani  parameters  d,k  as  computed  just  above.  The  exact  attri¬ 
tion  in  such  a  process  is  given  by 

AR  *  R^l  -  (1  -  £  [1  -  (1  -  d)R])B)  . 

Using  the  same  attrition  apportionment  as  in  (17),  these  three 
equations  yield  the  following  values  for  the  expected  attrition  to 
type  j  targets 

(18)  ARj  =  R  Vl  -  (1  -  |  [1  -  (1  -  d)R])3)  . 


c.  Prior  Apportionment  of  Searchers.  We  might  also  model  the 
attrition  process  as  a  number  of  smaller  engagements  by,  prior  to 
attrition  assessment,  dividing  the  searchers  among  different  types 
of  targets  on  the  basis  of  relative  numbers  of  targets  present  and 
vulnerable.  That  is,  type  j  targets  would  be  vulnerable  to  only 


B(i,j)  =  Bi 


type  i  searchers, 
interpret  this  as 
of  target. 


rather  than  to  all  B..  type  i  searchers.  One  can 
assigning  to  each  searcher  one  and  only  one  type 
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■'T7- 


T“. 


Using  (15),  we  see  that  attrition  to  type  j  targets  is  given  in 
this  case  by 


(19) 


d.  .) 

ii ' 


^B(M 


But  there  are  other  methods  of  prior  searcher  allocation  which 
are  more  closely  related  to  our  model  and  to  methods  developed  by 
L.  B.  Anderson  (1973).  The  point  is  to  compute  explicitly  the 
searcher  allocation  for  a  "typical”  target  set  (we  leave  the  term 
purposely  vague)  and  to  use  this  allocation  to  allocate  searchers 
in  each  attrition  computation,  as  described  below. 


Suppose  that  t.  is  the  proportion  of  a  "typical”  target  set 
3  t 

which  are  type  j  targets  and  that  a^  is  the  fraction  of  the  shots 
fired  by  a  type  i  searcher  which  are  directed  at  type  j  targets 
when  the  target  set  is  "typical”.  If  we  put 


(20) 


then  the  fraction  of  shots  fired  by  a  type  i  searcher  which  are 
aimed  at  type  j  targets  when  the  target  set  consists  of  targets 
of  type  j  (j  =  1,  ...,  N)  is  given  by 


-  -ikRk 


y 


and  we  could  then  define 


B(i, j ) 


a-  .R. 
-  r  3 

Bi  e  a 


A 


y 


to  be  used  in  (19). 
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An  interesting  question  at  this  point  is,  Can  the  a|j\  be  derived 
from  other  given  quantities  using  the  assumptions  set  forth  above? 

To  begin,  we  note  that  the  object  is  to  derive  a  distribution  of  the 
fire  of  each  searcher.  For  the  model  described  by  Al)  -  A7)  we  have 
the  following  result. 


(21)  PROPOSITION.  For  each  i  and  j  ,  the  conditional  probabi¬ 
lity  that  a  fixed  type  i  searcher  attacks  some  type  j  target,  given 
that  the  searcher  makes  an  attack,  is 


(22) 


1  - 


R.d. ,p(i, j ) 

3  2-3  ’ 

N 


n 

k=l 


(1  -  d 


ik 


> 


/ 


with  p(i,j)  defined  by  (11). 

PROOF.  In  the  proof  of  (9)  we  showed  that  the  probability  that  the 
kth  searcher  of  type  i  attacks  the  j£th  target  of  type  j  is  given  by 

p{Ai:j(k,A)}  =  di:jp(i,j)  . 


The  probability  that  the  k^  type  i  searcher  makes  an  attack  is 
one  minus  the  probability  that  he  fails  to  make  a  detection;  by  AS) 
this  latter  probability  is 


N 

n 

k=l 


(1 


R, 


d.,  ) 
ik' 


By  elementary  conditional  probability  arguments  (22)  is  now 
obtained.  IF 

Suppose  now  that  a  "typical"  target  set  consists  of  R^  targets 
of  type  j,  j  =  1,  ...,  N.  We  can  then  define 


t 


j 


»  t 
S  \ 
k=l  K 


> 
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and  uefine,  through  (22), 


.  R^d. •pt(i, j ) 

t  _  3  3Q 

ij  "  N  Rt 

i  -  n  a  -  d  A 

k=l 

t  t 

with  p  (i,j)  computed  by  (11)  with  the  there  replaced  by  R^,  Then, 

to  compute  attrition  i n  an  engagement  with  arbitrary  numbers  Rp  . . . , 

R^  of  targets,  we  would  compute  the  aik  by  (20)  and  apply  equation 

(19).  In  use,  this  procedure  would  require  only  one  computation  of 

P(i)j)’s — that  for  the  "typical”  target  set. 

The  rationale  behind  this  procedure  is  that  in  a  statistical  sense, 
of  the  long  run,  the  probability  that  a  given  type  of  target  is  fired 
upon,  given  that  there  is  a  shot  fired,  is  the  same  as  the  fraccion  of 
shots  fired  at  type  j  targets  if  the  autrjtion  process  were  carried 
out  experimentally  a  large  number  of  times. 

5.  A  Multiple  Shot  Model 

The  preceding  discussion  shows  that  the  complexity  of  (10)  and 
(11)  results  almost  entirely  from  the  Assumption  A3)  that  each 
searcher  fires  at  most  one  shot.  If  we  replace  A3)  by 

A3)'  Each  searcher  has  sufficient  firing  capability  to  fire 
once  and  only  once  at  each  target  he  detects, 

then  we  have  the  following  result. 

(23)  PROPOSITION.  Under  Assumptions  Al),  A2),  A3)',  A4)  -  A7), 

/  M  Bi\ 

(24)  AR.  =  E.  ^  (1  -  di.ki:j)  ) 

for  each  j  . 

We  omit  the  simple  proof.  The  equation  (24)  is  a  reasonable 
and  more  easily  computable  alternative  to  (10),  especially  for 
modeling  situations  in  which  the  one-shot  hypothesis  A3)  seems 
unrealistic.  oc 


6.  Exponential  Approximations 

As  in  the  homogeneous  case  various  exponential  approximations 
to  the  attrition  equations  of  this  section  are  possible,  some  of 
which  we  give  below.  Corresponding  to  (16)  are 

AEj  ~  Rj[l  -  exp  (-  i  Bik..[l  -  (1  - 


and 


t/  M  P  -d .  .R . - 

1 '  exp  I'  f,-  £  VijL1  -  e  ’]). 


Approximations  of  (17)  are 


and 


*3  ~  Ri[x  -  «P  (-  5  j,  W1  -  (1  -  di)R])] 

AR.~R.[l-SXp(-|  j^^l-  E))J  • 


Equation  (18)  is  approximated  as 
Ar 


-R^  1  -  exp(-  Si  [1  -  (1  -  d)R]j 


and 


AR.  ~  R 
3  3 


n  ^  /  Bk  -dR,\ 

.  1  -  exp(  -  [1  -  e  ]  J 


Finally,  two  approximations  to  (19)  are 


and 


iRj  ~  Rj[l  -  e*p(  -  i  B^l  -  U  -  dy)]R3  J 

AR.  ~  Rj[l  -  |  B^l  -  e  ^  ))]  . 
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Such  relations  as  these  should  be  used  with  great  care,  if  at 
all,  since  they  are  three-  and  four-tier  approximations  to  the 
correct  computation  given  by  (10)  and  (11). 
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IV.  AREA  FIRE  MODEL 


I.  The  Basic  Model 

We  give  in  this  section  an  axiomatic  treatment  of  a  classical 
area  fire  problem,  following  the  format  of  the  preceding  sections. 

Karr-  (1972b)  is  the  basis  for  this  presentation.  In  our  first  model 
here  both  defenders  and  incoming  shots  are  uniformly  distributed  over 
a  disk  in  the  plane;  it  is  shown  that  under  carefully  stated  assumptions 
an  independence  assertion  usually  thought  to  be  true  fails  and  we 
denote  some  attention  to  the  consequences.  The  second  model  is  an 
elementary  Markov  process  analysis  of  a  situation  in  which  the 
defenders  regroup  after  every  shot  so  as  to  defend  an  area  propor¬ 
tional  to  their  numerical  strength,  the  constant  of  proportionality 
being  fixed. 

2 

For  each  r  >  0  let  S  be  the  disk  in  R  with  center  0  and 
—  r  ~ 

radius  r  ;  we  denote  the  unit  disk  by  A.  Defenders  (or 
targets)  located  in  A  are  attacked  according  to  the  following 
assumptions . 

Al)  Each  incoming  shot  has  a  circular  area  of  (theoretical  -  see 
A3)  below)  lethality  whose  radius  is  uniformly  distributed  on  the 
interval  Tr^,  r^],  where  r 2  are  numbers  with 

0  <  rx  <  r2  <  1  , 

with  the  understanding  that  if  r^  =  r 2  the  radius  is  ^  almost  surely. 

A2)  The  center  of  the  area  of  lethality  of  each  shot  is  uniformly 

distributed  on  S.,  and  is  independent  of  the  radius  of  lethality, 
l-rg 

This  means  that  the  attacker  is  able  to  (and  does)  aim  his  shots 
to  the  extent  that  no  lethal  area  falls  outside  A  ,  but  that  centers 
are  otherwise  uniformly  distributed.  Note  also  that  no  point  of  A  is 
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Preceding  page  blank 


a  priori  immune,  although  points  near  the  boundary  of  A  appear 
safer. 

A3)  There  is  a  measurable  function  g  from  A  x  S,  x  TO,  r0] 

l~r2  * 

into  TO,  1]  with  the  interpretation  that  g(x,u,r)  is  the  prob¬ 
ability  that  a  defender  at  x  e  A  is  killed  by  a  shot  with  lethal 

radius  r  centered  about  u  e  S,  . 

l~Z2 

A4)  Each  defender’s  position  is  uniformly  distributed  on  A  and 
is  independent  of  all  incoming  shots. 

The  purpose  of  the  function  g  is  to  take  account  of  the  effect 
of  terrain,  shelters,  etc.,  in  decreasing  the  lethality  of  a  shot. 
Thus  one  may  suppose  that  the  lethality  of  Al)  is  hypothetical  in 
the  sense  of  representing  lethality  on  a  treeless  plain,  so  that 
g  has  the  effect  of  transforming  a  "potential”  lethality  which  is 
a  property  only  of  the  weapon  into  an  "actual”  lethality  which 
depends  on  both  the  weapon  and  the  physical  situation. 

Some  examples  are  helpful  at  this  point. 

EXAMPLES.  1)  the  function 

1  for  |x-u|  <  r 

(25)  g(x,u,r)  = 

0  otherwise, 

2 

where  |  •  |  denotes  Euclidean  distance  on  R  ,  is  the  classical 
”cookie  cutter”  damage  function:  anything  in  the  lethal  range  of 
the  shot  is  annihilated  with  probability  one  and  everything  not  in 
the  lethal  range  is  unaffected. 

2)  Suppose  the  subset  T  of  A  is  sheltered  in  the  sense  that 
no  defender  in  the  shelter  can  be  killed  except  by  a  direct  hit  on 
the  shelter.  We  would  then  have  for  x  e  T 

1  if  u  e  T  and  |x-u|  <  r 

g(x,u,r)  = 

0  otherwise  , 

while  for  x  j:  T,  g(x)  could  be  defined  by  (1)  or  by  some  other 
recipe. 
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Note  that  A4)  does  not  require  that  positions  of  different 
defenders  be  independent  of  one  another,  but  only  that  each  position 
be  uniformly  distributed  on  A  .  The  independence  assumption  will 
be  imposed  in  the  second  model. 

Let  us  consider  the  probability  p(l)  that  a  single  defender 

with  position  X  survives  a  single  shot  with  (center,  radius) 

pair  (C,R).  Writing  S  for  S,  ,  we  have 

i-r2 

(26)  1  -  p(l)  =  E  T g(X,C,R)]  = 

1  *2 

=  [  |  A j  •  |  S  j  •  (r2  -  r1))  1J'  dx  J  du  J  dr  g(x,u,r)  . 

A  S  r 

2 

where  I F I  is  the  two-dimensional  Lebesgue  measure  (=  area)  of  F  c  R 

o 

Of  course,  |a|  =  tt  and  |s|  =  n(l~  rg)  ,  but  we  leave  these 
quantities  unevaluated  for  the  sake  of  easier  interpretation.  It 
should  be  noted  that  throughout  this  section  u  and  x  are  two- 
dimensional  variables,  while  r  is  one-dimensional. 

(27)  PROPOSITION.  If  R  is  any  nonnegative  random  variable 
bounded  above  by  r2,  if 

p|c  e  I  ^  ’ 

2 

and  if  g  is  of  the  form  (25),  then 

(28)  i  -  p(i)  =  £IilL1  . 

|a| 

PROOF.  Let  F  be  the  distribution  of  R  and  G  that  of 

C;  then  with  the  assumption  on  g  ,  and  an  application  of  Fubini’s 
Theorem,  we  have 

r2 

1  -  p(l)  =  -^r  /  F<dr)  /  G(du)  J  dx  Ir(J  ^( | x-u| ) 


(here  I  is  the  indicator  function  of  the  set  B:  In(y)  =  1  if 
B  B 

y  e  B  and  0  otherwise) 


J  F(dr)  j  G(du)  ||x  e  A 
0  S 


But,  since  for  u  e  S  and  r<u,  s  (u)cA  , 

[  jx  e  A  :  |x  -  u|  <  r|  |  =  n  r2 
so  that  the  last  expression  is  equal  to 


1 


A 


/  F(dr) 
0 


2 

tt  r 


/  G(du)  = 
S 


Efn  R2] 

|A| 


5 


When  R  has  the  distribution  specified  in  Al), 


E[tt  R2] 


dr 


+  rir2 


+ 


As  a  special  case  of  (27)  we  may  take 

G({ 0})  =  1  , 


and  notice  that  the  probability  the  defender  survives  a  single  shot 
whose  center  is  fixed  at  the  origin  is  the  same  as  the  probability 
that  he  survives  a  shot  with  the  same  radius  of  lethality  distri¬ 
bution  whose  center  is  uniformly  distributed  on  S  . 

Let  p(n)  be  the  probability  that  a  single  defender  survives 
n  independent  and  identically  distributed  shots.  Then  it  is 
not  true  that 
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(29) 


P(2)  =  pCD2 


since  the  events 


and 


=  j defender  survives  first  shot} 


Eg  =  |  defender  survives  second  shot} 


are  not  independent,  but  only  conditionally  independent  given  the 
defender^  position  X  . 

Thus 


p( 2)  =  p|e1  0  Eg} 

=  Erp^  n  Egjx}i 
=  E[P|E1|X}  .  P|Eg|x}]  . 

It  follows  from  (26)  and  (27)  that  for  i  =  1,  2 

P{E.|X|  =  h(X) 

where 

-1  r  r2 

h(x)  =  (|s|  •  ( Tg  -  rx))  1  J  du  J  dr  I(r>c#)(|x  -  u|) 


1 

S| 


|u  e  S  :  |x  -  u|  >  r2}  | 


/ 

r]<l x_ui<r2 

ueS 


r2  '  rl 


( | x  -  u |  -  ^ )du 


Note  that  h(x)  is  the  probability  that  a  defender  at  x  s  A  survives 
a  shot  which  is  uniformly  distributed  on  S  .  It  follows  that 

(30  p(2)  =  |A|"1/h(x)2  dx 

n 
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and,  by  induction,  that 


(3.1)  p(n)  =  | A | -1  f  h(*)n  dx  . 

A 

The  probability  q(n)  that  a  single  defender  fails  to  survive  n 
independent  and  identically  distributed  shots  is  given  by 

q(n)  =  1  -  p(n) 

=  1  -  lAf1/ h(x)n  dx 
A 

=  |A|_1y* ri  -  h  (X)n]  dx  , 

A 

which  is  not  equal,  in  general,  to 

1  -  (1  -  q(l))n  =  1  -  (Iftf1/ h(x)  dx)n  . 

A 

The  expected  attrition  to  in  defenders  whose  distributions  satisfy 
A4)  is  then  given  by 

(32)  m  .  q(n)  =  mjAl"1^  (1  -  h(x)n)  dx  h  A(m,  n)  . 

A 

Let  us  again  point  out  that  A4)  is  satisfied  and  hence  (32)  is  valid, 
without  assuming  that  defenders’  positions  be  independent  of  one 
another.  For  example,  if  a  single  defender  chooses  a  uniformly  dis¬ 
tributed  position  and  other  defenders  position  themselves  at  the  same 
distance  from  0  in  such  a  way  as  to  make  all  angular  distances 
between  adjacent  defenders  equal,  then  A4)  is  satisfied.  Or  all 
defenders  may  be  placed  at  a  single  point  which  is  uniformly  distri¬ 
buted  on  A  ,  so  a  defender  seeking  to  minimize  his  expected  attrition 
gains  nothing  by  dispersing  his  personnel,  in  a  situation  modeled  by 
Al)  -  A4).  The  same  would  not  be  true,  however,  of  a  defender  who 
tries  to  minimize  the  probability  of  annihilation  of  all  his  personnel 
(provided  his  forces  are  more  than  one). 


2.  Edge  Effects 


As  will  be  discussed  iri  considerable  detail  below,  the  difficulties 
in  the  failure  of  (29)  arise  in  the  fact  that  h  is  not  constant  over 
all  of  A  :  clearly  h(0)  <  1,  while  h(x)  =  1  for  x  in  the  boundary 
Fr(A)  of  A  .  The  reason  h  is  not  constant  is  because  shots  are  aimed 
so  that  no  lethal  area  ever  falls  outside  A  ,  making  it  less  likely 
that  a  shot  whose  center  is  uniformly  distributed  on  S  will  kill  a 
defender  whose  position  is  close  to  Fr(A).  This  is  an  edge  effect. 


It  will  be  proved  in  (38)  that  h  is  constant  over  most  of  A  ; 

we  will  then  consider  the  effect  of  neglecting  the  set  where  h  ^  h(0). 

Another  alternative  which  makes  h  constant  over  all  of  A  would  be 

to  make  shot  centers  uniformly  distributed  on  S,  ,  but  if  this  is 

l+rj’ 

done  then  (28)  is  no  longer  true.  The  difficulties  arise  in  "edge 
effects"  near  rhe  boundary  of  A  .  One  is  tempted  to  conclude  that 
edge  effects  can  be  neglected  if  r2  is  small;  tc  a  certain  extent  this 
is  so,  or  that  these  edge  effects  are  a  result  of  the  formulation  of 
the  model  rather  than  the  underlying  physical  process.  We  believe  the 
second  conclusion  is  unjustified.  In  physical  models,  from  fluid  dyna¬ 
mics  to-  statistical  mechanics,  of  finite  systems  it  is  precisely  the 
boundary  behaviors  which  must  be  considered  most  carefully  and  which 
make  the  theories  valuable;  we  feel  the  same  is  true  here. 


3.  Repositioning  in  the  Same  Area 

We  now  proceed,  as  promised  above,  to  a  detailed  analysis  of  the 
failure  of  (29)  and  one  way  of  rectifying  this  difficulty  (if  one 
considers  it  to  be  a  difficulty)  by  imposing  an  additional  (but  not 
very  plausible)  assumption. 

Suppose  we  add  to  Al)  -  A4)  the  following  hypothesis: 

A5)  Shots  are  made  one  at  a  time  and  between  each  pair  of  shots 
every  defender  repositions  himself  according  to  a  uniform  distribution 
on  A  ,  independent  of  the  past  positions  of  all  defenders  and  of  the 
entire  process  of  incoming  shots. 
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I 

« 


— 


We  begin  by  considering  the  survival  probabilities  for  a  single 

defender;  suppose  there  are  n  shots  that  his  positions  (provided 

he  lives  long  enough  to  occupy  them)  are  X,,  , X  and  that  (C. ,  R  ) 

X'  n  x  1  ’ 

•••»  (Cn,  ^n)  are  (center,  radius)  pairs  of  the  incoming  shots. 
Continue  to  assume  g  is  of  the  form  (25).  Then  for  k  =  1,  ...,  n 
the  probability  of  surviving  the  first  k  shots  is 

p<k>=pjjUxs  *  vCj)H 

Hence  under  A5)  and  the  assumptions  of  Proposition  (27) 


(33) 


P( 


n)  =  /l  -  Sins!!)" 

\  |A|  / 


n  =  1,  2,  ..., 


where  R  is  a  random  variable  with  the  same  distribution  as  R,,  . . . ,  R  . 

1’  n 

From  (33)  it  is  evident  that 


(34) 


q(n)  =  1 


2l\ 


(l  ETttRII  . 

\  |A|  / 


=  1  -  (1  -  q(l))n. 


We  emphasize  that  (34)  is  valid  only  under  the  restrictive  and  rather 
implausible  hypothesis  A5)  and  that  A5)  is  meaningless  in  a  situation 
where  r.  shots  are  fired  simultaneously  and  independently  from  differ¬ 
ent  weapons. 

For  small  values  of  r ^  and  large  values  of  n  ,  one  may  make  the 
approximation 


(3f) 


p(n)  ~  exp  (-  nE[rrR2]/|A| ) 


although  neither  its  desirability  nor  its  necessity  is  clear,  since, 
after  all,  exact  expressions  such  a?  (33)  are  available.  If  (35)  is 
employed,  the  expected  attrition  A(m,  n)  to  m  defenders  by  n  shots 
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r. 


is  approximately 


A(m,n)  ~  m(l  -  exp  (  -  nE{>R2]/| A|  ))  ; 


the  true  value  is 

(36)  A(m,  n) 


It  seems  reasonable  to  define  the  quantity  nEruR^'J  as  the  (expected) 
"antipersonnel  potential"  of  n  shots.  It  is  the  word  "potential" 
which  is  crucial:  this  is  the  expected  lethal  area  of  n  shots  with 
the  same  radius  of  lethality  distribution  whose  areas  of  lethality 
do  not  intersect  (which  is  the  best  situation  the  attacker  can  hope  for). 


Continuing  to  assume  Al)  -  A5)  and  the  hypothesis  of  (27),  we  see 

that  the  probability  p(n1}  . ..,  n^)  that  a  single  defender  survives  n.^ 

independent  shots  from  a  type  1  weapon  with  expected  area  of  lethality 
2 

E[ttR^]}  n„  independent  shots  from  a  type  2  weapon  that  are  independent 

of  the  shots  from  the  type  1  weapon  and  have  expected  area  of  lethality 
o 

ECttR^L  and  so  on,  is 

n. 

l 

j  /  E[ttR?1 

(37)  p(n  ...,  n  )  =  n  1  -  — - — 

3  i=l\  |A| 

We  may,  provided  n^,  n^  all  be  large,  approximate  the  corresponding 

kill  probability  q  =  1  -  p  by 


q(n1,  ru)  =  1  -  pCn1}  ...,  n.. ) 

~  1  "  e-P  —  s  niE^R?])  • 

\  | A|  i=l  1  1  / 

2 

Again,  it  may  be  reasonable  to  interpret  E  r^EfttR^]  as  antipersonnel 
potential.  It  should  be  noted,  however,  that  no  such  scalar  quantity 
appears  in  the  exact  equations  (34)  and  (37),  let  alone  the  equations 
derived  without  the  use  of  A5). 
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4.  The  Effect  of  Edge  Effects 


An  alternative  analysis  of  the  failure  of  (29)  is  based  on  the 
following  result.  We  assume  until  further  notice  that  r2  <  1/2. 

(38)  PROPOSITION.  If  |x|  <  1  -  2r2,  then 
h(x)  =  h(0) 

~  “ — “ — ^2  [t1  "  2r2)  +  2/3  (r2  +  rir2  +  rl'  '  ri(r2  +  rl^] 

PROOF.  For  |x|  <  1  -  2r2  we  have,  from  the  computation  of  h  given 
previously, 


h(x)  =  -—-[({ u  e  S  :  |x  -  u|  >  r2}| 


/  (|x  -  u|  -  rp  du] 

1  r!<l x“ul<r2 


ueS 


2. 

TTr0) 


-nr  /  (ix  ■  ul  _  ri)  du]  * 

1  r1<!x-u|r2 
|u|<1-r2 


In  the  second  term,  make  the  change  of  variable  y=u-  x  (y  is  two- 
dimensional)  to  obtain 


- - -  f  (|y|  “  *x)dy  • 

! S I ( r2  "  rl^  r^<jyj<r2 
|x+y|<l-r2 
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For  | yj  <  r2  and  jx|  <  1  -  2r2, 

|x  +  y|  <  |xj  +  |y|  <  1  -  r2  , 

so  the  second  restriction  on  the  domain  of  integration  is  irrelevant. 
Hence  for  |x|  <  1  -  2r2  , 


h(x)  =  h(0) 


1  f 
~  Is!  L 


2  2 
tt(1  -  r2)  -  nr2 


■T-rT  f  (hi  -  V  dy] 

ri<|y|<r2 


— ?  I*1  -  2V  +  KT^Fr  /  (Wi-V 

-  r„j  L  •  z  1  J 


dy 


tr(l  -  r2) 


We  may  evaluate  the  integral  changing  to  polar  coordinates  y  -  (r  cos  9, 

-  1 

r  sin  0),  where  r  =  |y{  and  Q  =  tan  (yg/y^),  with  the  result 


f  (|y|  “  t1)  dy 

ri<|y|<r2 


2tt 


-L-—  J  J  (r  -  rx)  r  dr  dg 


r1<r<r2  0 


2tt 


r2  -  r 


-  J  (r2  -  r^r)  dr 


r^<r<  r2 


_  2n 


hr  [V3(r|  -  4)  -  ^(r22  -  ^)] 


which  completes  the  proof. 
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COROLLARY. 


1 1 x  :  h(x)  ±  h(0)}j  =  4rr(r2  -  r^)  . 

PROOF.  This  requires  only  noting  that  h(x)  is  strictly  greater  than 
h(C)  if  | x J  >  1  -  2r2  and  a  simple  computation.  Q 

Since 

n 


(1  -  q(l)) 1  = 


n=fe-(hcx)dx) 


while 

q(n)  =  1  -  —  f  h (x)n  dx  , 

lSl  { 


we  see  that  the  failure  of  q(n)  to  equal  1  -  (1  -  q(l))n  under  Al)  - 
A4)  and  the  hypotheses  of  (27)  is  due  to  the  inequality 


—  f  h(x)ndx  ? 

m 

We  know  that  if  G  is  a  probability  measure  and  f  =  c  is  a 
constant  function  then 

J  fn  dG  =  ( Jf  dG)n  =  cn 


for  all  n  >  0. 


Since  we  showed  in  Proposition  (38)  that  h  is  constant  except 
over  a  subset  of  A  of  small  Legesgue  measure  (provided  r„  is  small) 


it  seems  plausible  that 


to 


/>• 


should  serve  as  an  approximation 


This,  in  fact,  is  so  in  two  senses,  which  we  note  next. 


\-C 


— 


% 


(39) 

PROPOSITION. 

For  any 

fixed 

n 

lim 

1 

■  f  h  (x)n 

dx 

1 

r2-0 

.  iA! 

J 

A 

\ 

1*1 

and  for 

any  fixed 

r2 

>  o, 

‘ 

lim 

1 

-  f  h(x)n 

dx  -( 

1 

n-*<= 

A 

J 

A 

\ 

1 A  | 

=  0  , 


rt 


=  0  . 


We  omit  the  proof. 

COROLLARY .  For  each  fixed  n  >  1  and  each  6  e  TO,  1), 
JL  f  h(x)n  dx  =  h(0)n  +  o(4)  ,  r  -  0 

wi 


and 


n 

/_i_  j h(x)  dx  j  =  h(0)n  +  o'.r|)  , 

\|AI  a  / 


r2  -  0  ’ 


The  qualitative  content  of  the  corollary  is  that  for  fixed  n^  and 
small  r2  the  error  in  approximating  either  q(n)  cr  1  -  (1  Q(l))  by 

1  -  h(0)n,  and  hence  the  error  in  the  approximation 


n 


q(n)  ~  1  -  (1  -  q(l» 

is  (roughly)  proportional  to  r2;  estimates  of  the  constant  of  propor¬ 
tionality  may  be  made  if  required. 

Let  f(x)  =  1  -  h(x),  so  that  f(x)  is  the  probability  that  a  defender 
at  x  e  A  is  killed  by  a  single  shot  whose  center  is  uniformly  distri¬ 
buted  on  S,  and  is  independent  of  the  radius  of  lethality  which  is 
l-r2 

uniformly  distributed  on  [r^  r2].  Then  the  approximation 
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q(n)  ~  1  -  h(0)n 

may  be  written  as 

q(n)  ~  1  -  (1  -  f(0))n  . 

Since  both  sides  here  go  to  zero  as  n  -  ®,  we  may  make  the  further 
approximation 

q(n)  ~  1  -  exp  (  -  nf(0)) 

which  is  valid  for  small  r^  and  large  n  . 

5.  Repositioning  With  Fixed  Density 

We  will  next  discuss  a  generalization  of  the  model  based  on  Al)  - 
A4)  which  treats  incoming  shots  as  occurring  one  after  another  and 
reqjires  that  the  defender  reposition  his  forces  after  each  shot  in 
such  a  way  as  to  maintain  a  fixe i  density  of  forces.  Specifically, 
we  will  assume  Al)  -  A4),  that  there  are  initially  m  defenders,  and 
the  following  additional  hypotheses: 

Ab)  Let  p  =  m/|A| .  If  after  a  given  shot  there  are  k(l  <  k  <  m) 

surviving  defenders,  each  of  those  k  defenders  positions  himself  on 

the  disk  S  with  radius  r(k)  =  (k"pn)2  =  (k/m)2  and  center  0  , 
r(k)  . 

independent  of  the  repositionings  of  the  (k  -  1)  other  survivors,  o~  che 

past  positions  of  all  defenders  and  of  the  entire  process  of  incoming 
shots . 

A7)  Initial  positions  of  the  m  defenders  are  independent  of  one 
another  and  of  all  incoming  shots . 

Hence  the  density  of  the  repositioned  defenders  (over  the  smaller 
set  S  »  to  which  they  have  retreated)  is  k/nr(k)  =  p  so  the  density 
of  defenders  (relative  to  the  set  over  which  they  are  uniformly  dis¬ 
tributed-  -not  to  A  !)  remains  constant. 

For  simplicity  we  further  assume  that  g  is  of  the  form  (2_>), 
that  all  shots  have  lethality  radius  r  almost  surely,  that  the 
attacker  is  unaware  of  the  pullback  of  defender  personnel  and 
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continues  to  fire  shots  whose  centers  are  uniformly  distributed  on 
S^_r>  ar*d  that  r  <  1/m.  Let  Tn  be  the  number  of  defenders  remaining 
alive  after  the  n^  shot:  The  important  result  concerning  the  sto¬ 
chastic  process  (T  )  x  is  the  following. 

(40)  THEOREM.  Under  Al)  -  A4),  A6),  A7),  and  the  additional 
assumptions  stated  in  the  preceding  paragraphs,  T  =  (Tn)n>1  is  a 
Markov  process  with  state  space  E  =  |  0,  1,  . ..,  m},  initial  distri¬ 
bution  a  given  by 

(41)  aOO  =  P\'S1  =  k|  =  1  -  r2)k(r2)m  k  ,  k  =  0,  ...,  m 

and  transition  matrix  P  given  by 


(42a)  F(0,  j)  =  6(0,  j)  ,  j  =  0,  1,  ...,  m 

while  for  1  <  i  <  m 


(42b) 


P(i,  j) 


f.(u))j  f^u)1^  du  , 
if  D  b,  ...,  i 

if  j  =  i  +  1,  . . . ,  m  , 


where  S  f(0,  u)]  is  the  disk  with  radius  r  and  center  (0,  u)  and 


_  !Sr(i)  n  u»i 

;vUJ  "  O  > 


TTV(i) 


and  finally 


(42c)  P(m,  j)  =  a(j)  ,  j  =  0,  ...,  m  . 

PROOF.  As  in  the  derivation  of  (31)  one  must  guard  against  unwarranted 
independence  assumptions.  That  ( 42a ) holds  is  clear  from  physical 
interpretation. 

If  X1?  . . . ,  X  arc  the  positions  of  the  defenders  prior  to  the 
first  shot  and  the  position  or  the  center  of  the  first  shot  then 


P{T1  =  k|  k  (.^)p|xi  i  S^),  ...,  S^), 

Xk+1  e  Sr^Cl^  Xm  e  Sr^Cl^ 

(because  defenders  are  indistinguishable) 

=  (£)  E[P{X1  i  Sr(C1),  ...,  Xk  e  S^), 

Xk+1  e  W’  Xm  e 
C  ro-k 

='(”)  ~r  J  ^  ?{h  t  VW  P1\+l  e  SM 


s 


(because  X,,  . ..,  Xm  are  independent  of  one  another  and  of  C.^) 


k  m-k 

■(; K1-2^)^) 


proving  (41)  and  (42c).  A  similarly  pleasant  simplification  in  the 

proof  of  (42b) is  not  possible  because  a  shot  uniformly  distributed 

on  S  does  not  have  a  probability  of  killing  a  defender  uniformly 
1-r 

distributed  on  S  ,.N,  which,  given  the  center  of  the  shot,  xs  con- 

r(i) 

stant,  as  was  true  above,  except  for  i  -  m. 

Suppose  C,  C2,  C^+1  are  the  centers  of  the  first  k+1  shots 

and  let  X, ,  . . . ,  X™  be  the  positions  assumed  by  the  Tk  defenders 

1  k  .  .  ! 

surviving  the  first  k  shots.  Since  by  A6)  X,,  X™  are 

k 

independent  and  identically  uniformly  distributed  on  Sr^  given 

Tk,  and  are  independent  of  C^,  . . . ,  C^,  and  of  all  previous  defender 

positions,  it  follows  that  XT  are  conditionally  independent 

k 

of  T  ,  T  •  given  Tk,  and  hence  that  the  Markov  property  holds: 

piTk+l  =  ^  lTi*  •••»  =  P^Tk+l  " 
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It  remains  to  prove  that 


plTi»i  =  3lTkl=  p<Tk’  i) 

with  P  defined  by  (42);  we  may  and  do  assume  that  >  1.  Noting 

that  P{T^+1  =  j|Tk}  is  the  conditional  probability  given  that 

exactly  -  j  of  the  points  . ..,  X^  are  within  the  disk  Sr(Ck+1) 

k 

and  again  using  the  indistinguishability  of  the  defenders,  we  obtain 
P<Tk+l  =  ^Tk  = 

=  (/)  P^X]L’  ^  Sr(Ck+1^  X3+1>  XTk  G  M^+l^W 

"  (j  )  E!^X1’  Xj  ^  VCk+l^’  Xj+1»  Sr^Ck+l^lTk’  Ck+lMTk^ 

’  _  (  J  /  ?k"j 

=  (/)  TT'l  /  dy  PK  1  Sr<y)l  T1  £  Sl<y)l 

1'r  sl-r 

By  rotational  symmetry  of  a  uniform  distribution  on  a  disk, 
pjX-^  e  Sr(y)}  depends  only  on  |y|  and  since  |Ck+1|  is  distributed  on 
[0,  1  -  r]  with  density  (1  -  r)-22u  the  last  expression  becomes 


j  .  j 

(Ik)  J  2udu  P^XX  i  Sr[(0,  u)]}  PjX1  e  Sr[0,  u)]j 

Since  X]L  is  uniformly  distributed  on  we  have 

,  lSr(i)  n  Sr[(0’  U)]l 

P{X1  e  Srf(0,  u)]f=  f.(u)  = - ^ - j - 


Vj 


r(i)1 


which  completes  the  proof. 
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Further  remarks  concerning  the  functions  are  in  Appendix  A. 
The  complexity  of  the  expressions  for  the  elements  of  P  and  the 
exact  form  of  the  f^  as  given  by  (A. 4)  seem  to  preclude  efficient 
computation  of  P  on  a  computer. 

Appendix  B  considers  this  process  with  an  exogeneous  time  scale, 
rather  than  that  determined  by  numbers  of  shots. 
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APPENDIX  A 


Let  us  consider  in  more  detail  the  functions  f.  defined  in  (40). 

To  find  the  exact  form  of  these  functions  one  must  compute 
!®r(i)  n  Sr[(0,  u)]|  for  u  e  [0,  1  -  r].  It  is  evident  that 

2 

nr  ,  if  0  <  u  <  r(i)  -  r 

|sr(ij  n  sr[<0,  U)]|  = 

0  ,  if  u  >  r(i)  +  r  , 

, 

so  that  only  u  e  (r(i)  -  r,  r(i)  +r)  remain  to  be  considered.  If  one 
draws  a  picture  with  the  upper  half  of  the  circle  bounding  which 

has  the  equation 

(A.l)  y  =  (r(i)2  -  x2)1/2  , 

and  the  circle  which  is  the  boundary  of  Srr(0,  u)]  and  whose  equation 
is 

2  2  2 

(A. 2)  x  +  (y  -  u)  =  r  , 

he  sees  that  conventional  methods  of  analytic  geometry  may  be  used  to 
evaluate  the  area  of  js^.^  H  Sr[(0,  u)]|. 

Solving  (A.l)  and  (A. 2)  simultaneously  shows  that  the  two  curves 
intersect  at  the  points  (m^u),  m2(u)),  (  -  m^u),  m2(u)),  where 
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Let 


% 

8-j^x)  =  (r(i)2  -  x2)  ,  -  r(i)  <  x  <  r(i) 

2  2  ^ 

32(u,  x)  =  u  +  (r  -  x)  ,  -  r  <  x  <  r 

2  2  ^ 

83(u,  x)  =  u  -  (r  -x)  ,  -  r  <  x  <  r 

Note  the  implicit  dependence  of  8-^  on  i. 

2  9  ?5 

From  (A. 3)  we  see  that  m^(u)=r  when  u=(r(i)  -r  )2  and- that  m1(u)<r 
for  all  other  r  in  [r(i)-r,  r(i)+r]. 


2  2  ^ 

Thus  for  (r(i)  -  r)  <  u  <  r(i)  +  r  , 


m^u) 


m1(u) 


Isr(i)  n  sr  C(°j  u)]|=  J  B^Cx)  dx  "  f  8-(u,  x)  dx 


-m^u) 


2  2 

while  for  r(i)  -  r  <  u  <  (r(i)  -r) 


I Sr( i )  n  Sr 


m^u) 


’  U^*  =  f  B2(u> 

-m^u) 


m1(u) 

x)  dx  -  J  B-^C x)  dx 

-m1(u) 


Using  tables  of  integrals  we  obtain  the  following  evaluation  of 
the  functions  f^  . 
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APPENDIX  B 


The  regrouping  model  suggested  here  is  capable  of  handling  random 
numbers  of  shots  between  repositionings  of  the  defenders.  We  now  give 
the  theory  and  computations  necessary  to  make  this  generalization. 

Consider  the  Markov  process  (T  )n>1  with  transition  matrix  P  . 

Let  U  be  the  number  of  shots  on  day  n  ;  we  assume 
n 

A9)  Numbers  of  shots  on  various  days  are  independent,  identically 
distributed,  bounded,  positive,  integer-valued  random  variables, 
independent  of  the  positions  and  lethalities  of  all  shots  and  of 
defender  positions. 

If  we  put 


Vn  =  Ui+  ...  +un  ,  n  -  1 

Then  S  =  T  is  the  number  of  surviving  defenders  at  the  end  of  the 

th  a™  Yn 
n  day. 

(B.l)  PROPOSITION.  Under  the  preceding  assumptions,  (Sn)n>1  is  a 
Markov  process  with  initial  distribution  p  given  by 

H(k)  =  P\S±  = 

v  *0  m  T-l 

(B.2)  =  E  cp(i)  S  (j,  k) 

i=l  j=0 

and  transition  matrix  Q  given  by 

(B.3)  Q(ij  i)  =  E  cp(m)P  (i,  j)  . 

m=l 
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Preceding  page  blank 


Here  cp  is  the  common  distribution  of  U^,  Ugs  ...  and  is  the  upper 
bound  on  the  number  of  shots  per  day. 

PROOF.  According  to  the  terminology  of  Feller  (1966),  (Sn)  is  sub¬ 
ordinated  to  (T  )  by  the  distribution  co  .  The  proof  required  here  _ 

is  a  straightforward  application  the  theory  derived  in  Feller  (1966).  — 

Likely  candidates  for  c p  would  be  Poisson  and  geometric 
distributions . 
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